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We provide a new algorithm that translates a unitary matrix into a quantum circuit ac- 
cording to the G = KAK theorem in Lie group theory. With our algorithm, any matrix 
decomposition corresponding to type-AIII KAK decompositions can be derived accord- 
ing to the given Cartan involution. Our algorithm contains, as its special cases, Cosine- 
Sine decomposition (CSD) and Khaneja-Glaser decomposition (KGD) in the sense that 
it derives the same quantum circuits as the ones obtained by them if we select suitable 
Cartan involutions and square root matrices. The selections of Cartan involutions for 
computing CSD and KGD will be shown explicitly. As an example, we show explicitly 
that our method can automatically reproduce the well-known efficient quantum circuit 
for the n-qubit quantum Fourier transform. 
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1 Introduction 

Decomposing a unitary matrix into an efficient sequence of elementary gates is a fundamental 
problem in designing quantum circuits. There are two types of decomposition: One is exact 
decomposition where an arbitrary unitary matrix is decomposed precisely into a sequence of 
elementary gates, such as arbitrary single-qubit rotations and the CNOT. The other involves 
approximate strategies by which an arbitrary unitary matrix is decomposed approximately 
into a sequence of a fixed set of elementary gates, as shown in Solovay-Kitaev theorem (cf. [1], 
Appendix 3). In this paper, we treat exact decomposition. 

Cosine-Sine decomposition (CSD), which is a well-known algorithm in numerical linear 
algebra, was the first algorithm utilized for this purpose [2]. CSD applies the well-known 
algorithm for computing generalized singular value decomposition (GSVD). In CSD, we first 
divide an input matrix g into four square matrices and then apply SVD to each matrix. Then, 
we have 

(flu Si 2 \ = (ui 0\/c -s\U 0\ 

\fl21 522/ \0 u 2 J \s c J \ V 2 J 

where c = diag [cos(Ci), cos(£2), • • • , cos(C2»-i)]) s = diag [sin(Ci), sm(( 2 ),•••, sin(C 2 n-i)], 
then <7n = uicvi, 512 = — U1SV2, 521 = U2SU1, and 1722 = u 2 cv 2 are SVD. Recursively per- 
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forming decomposition JQl leads to a quantum circuit. The CSD-based algorithms are easy 
to implement on a computer because algorithms for calculating GSVD are well-known, and 
software libraries including GSVD are available. Some CSD-based algorithms [3-6] have been 
investigated with the aim of improving Barenco's result that an arbitrary 2™ x 2™ unitary 
matrix is composed of 0(n 2 4 n ) elementary gates [7]. And improvement to 0(4 n ) elementary 
gates has been reported by Mottonen et al. [4] and by Shende et al [6] . 

On the other hand, Khaneja and Glaser provided another kind of decomposition [8], which 
was later named KGD. KGD lies within the framework of the G — KAK theorem (cf. [9], The- 
orem 8.6) in Lie group theory. This theorem shows that an element g E SU (2™) is decomposed 
into matrix products k±ak2 for some k\, &2 € exp(t) and a <E exp(f)). Here, su(2") = t © m 
is a Cartan decomposition in Lie algebra su(2"), t and m = i 1 - arc orthogonal vector spaces 
contained in su(2"), and f) is a maximal Abelian subalgebra (a Cartan subalgebra) contained 
in m (cf. [15], §VI.2). Matrices hi, a, and ki are not uniquely determined from g. They 
depend on the selections of the bases of t, m, and f); besides, they are not determined even if 
bases are selected. Khaneja and Glaser provided a particular selection of bases of {, m, and f) 
in Ref. [8] so that the selection matches an NMR system, and they proved that a time-optimal 
control on a two-qubit NMR quantum computer can be obtained from the decomposition [10]. 
Thus, KGD can be regarded as the G — KAK theorem on the particular bases. It should 
be noted that KGD does not give a unique translation of the input matrix into a quantum 
circuit. 

Bullock [11] showed that CSD can also be regarded in the framework of the G — KAK 
theorem; i.e., CSD uses the type-AIII KAK decomposition with the global Cartan decompo- 
sition 9 defined as 6(A) = <7\ z X<j\ z for A S SU(2 n ), where <7j Z denotes that the operation 
defined as the Pauli matrix <j z acts on the j-th qubit. He also introduced a method that trans- 
lates matrices U, S, and V in Q into fci S exp(fi), a S exp(f)), and fc 2 € exp(t), respectively, 
where £ and f) are the ones defined in KGD. Here, KGD corresponds to G — KAK decomposi- 
tion with the selection of defined as 0(A) = a nz Xa nz . We can thus produce a KGD-based 
quantum circuit by combining Bullock's translation and the CSD-based algorithms. 

We introduce a new algorithm that translates a 2" x 2™ unitary matrix into a quantum 
circuit according to the G = KAK theorem. The algorithm can derive any matrix decomposi- 
tion corresponding to type-AIII KAK decompositions for the given global Cartan involution 
O. The algorithm contains, as its special cases, both CSD and KGD in the sense that it 
derives the same quantum circuits as the ones calculated by them if we select suitable Car- 
tan involutions and square root matrices. Here, we select 0(A) as <J\ z Xa\ z for CSD and as 
a nz Xa nz for KGD, where A € SU(2 n ). The strategy utilized in our algorithm is related to 
those used in Refs. [12,13]. However, those strategies provided methods for computing type- 
All KAK decomposition; no translation between type-All decompositions and type-AIII 
decompositions was provided. Furthermore, the method utilized in Ref. [12] is different from 
ours in the square root matrix calculations, i.e., methods for calculating m from m 2 (where 
g = km is a global Cartan decomposition of the input matrix g). In the method proposed in 
Ref. [12], first, a square root matrix is calculated in Lie algebra level. And then it is translated 
into an element in Lie group level via exponential mapping. In contrast, with our method, a 
square root matrix is calculated directly at the Lie group level. 

Although our algorithm contains CSD and derives any matrix decomposition correspond- 
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ing to type-AIII KAK decompositions according to the given Cartan involution, the efficiency 
for calculating a decomposition is not sacrificed. The reason is as follows: Roughly speaking, 
to decompose g into klak2, the CSD-based algorithms apply SVD to four 2 n_1 x 2™ _1 matri- 
ces (gtx, g\2, 321, and <?22 in Eq.Q), while our algorithm applies eigenvalue decomposition to 
2™ x 2™ matrix. Therefore, the efficiencies for computing SVD on four 2™~ 1 x 2™ _1 matrices 
and for computing eigenvalue decomposition on one 2™ x 2™ matrix are the same. 

In addition, our algorithm might have an advantage over CSD when we would like to 
determine a class of quantum circuits for a given class of matrices. The reason is as follows: 
In CSD-based algorithms, it is difficult to formulate a class of matrices u\, U2, V\, and 112 
such that relation holds for a given class of input matrices. Actually, to reproduce the 
well-known QFT circuit by using CSD [2,14], Tucci changes the rows and columns of the 
QFT matrices beforehand and makes each submatrix hold a convenient form, which can be 
written by the (n — l)-qubit QFT. It would not be possible to describe the general form of 
the decomposition when the input matrix does not have a convenient form like QFT. On the 
other hand, our algorithm does not require such a preliminary change of rows and columns. 
All matrices appearing through our algorithm can be described using the input matrix <?, 
the given global Cartan involution 0, and the eigenvalues and eigenvectors of these matrix 
products. This will be shown explicitly as an example of the QFT decomposition (Section 4). 

The paper is organized as follows: In the following section, we cover some preliminaries 
about notations, the G = KAK theorem, and KGD. Section [3] presents our algorithm for 
computing the decomposition follows from the G — KAK theorem. Section 0] presents de- 
compositions of the n-qubit QFT using our algorithm and CSD-based algorithms. We show 
that we can produce the well-known QFT circuit by using these matrix decompositions. 

2 Preliminaries 

2. 1 Notations 

Let a x , a y , and <j z denote the Pauli matrices and I® s be a 2 s x 2 s identity matrix (J = 2 1 x 2 1 ). 
We use aj a to denote the Pauli matrix acting on the j-th qubit; 

(a = x,y, or z). Let [/cnot denote the standard CNOT gate, H denote a Hadamard gate, 
and R X (C) = ex P( — iC a x)- All these notations follow those in Ref. [1]. 

2.2 G= KAK theorem 

The G = KAK theorem for compact groups (cf. [9], Theorem 8.6) provides a framework for 
decomposing g £ SU(2 n ) into the following matrix products: 



Here, su(2") = t © m is a Cartan decomposition, where t and m = are orthogonal vector 
spaces, and [) is a Cartan subalgebra, that is, a maximal Abelian subalgebra contained in m. 

Let 6 denote the Cartan involution of its Lie algebra su(2 n ); i.e., (i) 6 2 = I® n (6 ^ I® n ) 
and (ii) 9 is an automorphism of the Lie algebra su(2 n ). And let the global Cartan involution 
(cf. [15], p. 362) of SU(2 n ) be 0. Then I and m have the following property: 



g = k x ak 2 , 



ki,k 2 £ exp(t),a £ exp(()) C exp(m). 



(2) 




if x £ 6 
if x £m 




if X £ exp(«) 
if X £ exp(m) 



(3) 
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6 and exp(t) m 



Fig. 1. Patterns of an element of t, exp(t) and m for a three-qubit system, where each square 
represents an element of an 8 X 8 matrix. The white elements are always zero, and the black 
elements take some value that depends on the input matrix. 



Three types of t-algebra, named AI, All, and AIII, arise for su(2 n ). Here, AI, All, and 
AIII correspond to i = so(2"), t = Sp(2 n ), and s[u(p)®u(q)] (p + q = 2"), respectively (cf. [9], 
p. 518). 

2.3 Khaneja-Glaser decomposition (KGD) 

Khaneja and Glaser provided a particular selection of bases of 6, m, and f) (cf. [8], Notation 3 
and 5) so that the selection matches an NMR system. Notice that we use f) instead of f)(n). 
Here, generators of m, and f), are denoted to be tensor products of the Pauli matrices; 

t = span {A ®a z /2,B®I,i<T nz /2 \ A, Besu(2 n - 1 )}, (4) 
m = span {A(g> a x /2,B ® a y /2,ia nx /2,ia ny /2 | A, B G su(2™ -1 )}. (5) 

Here, generators of { and m have a specific operation on the last qubit; i.e., a z or / for 
generators of 6 and a x or a y for generators of m. Thus, to determine t and m, the Cartan 
involution 6 and the global Cartan involution can be chosen as follows: 

6(x) = <T nz xa nz , 9(1) = <7 nz Xa nz . (6) 

Since 9(a z ) = a z , 0(1) = I, 6(<j x ) = —o x , 0(a y ) = —o yi we can check that the above I and m 
satisfy relation © when 9 and 8 are chosen as I©- 

For the number of qubits n > 3, 1 and m have specific patterns, as shown in Fig.^ because 
all generators defined in 10} and © have these patterns. Note that in contrast to an element 
of exp({!) taking the same pattern as an element of t, an element of exp(m) does not take the 
same pattern as m. This property enables us to apply the KAK decomposition recursively, 
as shown in Fig. 

3 Our algorithm 

3.1 Basic ideas 

We provide a new constructive algorithm that computes a decomposition based on the 
G = KAK theorem. Here, we choose a Cartan subalgebra f), which is different from the 
t) used in KGD. Since Cartan subalgebras are Abelian, they can translate each other by 
f) = Ad Tecxp ( ( )(f)). Here, T is fixed for given 1) and f). It should be noted that G — KAK de- 
composition for the fixed input g G SU(2 n ) is not unique. We provide the following theorem. 
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Theorem 1 Let g G SU(2 n ) be the input matrix. If g has a global Cartan decomposition 
g — km (k £ exp(t), to € exp(m)), then to 2 is uniquely determined by m 2 = Q(g <! )g. 

Proof. From @(g')g = Q(rrvk')km = Q(rrv) Q(k^)km = mk^km — to 2 . □ 

Theorem 1 shows that the fixed global Cartan involution only determines to 2 . There- 
fore, arbitrariness remains in the selection of to, and also k. Furthermore, a in J5J) has also 
arbitrariness because it follows from a decomposition of to, to = k^ak, where k € exp(£). 
(See, [15], §AII.3). Using Theorem 1, we compute the decomposition in J5J as follows: 

1. Compute to 2 = ®(g^)g- 

2. Decompose to 2 = pbp^ such that p € exp(l) and b £ exp(t)). 

Such decomposition always exists because m 2 € exp(m) (cf. [15], Proposition 7.29). We 
compute p using eigenvectors of to 2 . We show examples in Section 13.21 

3. Find y such that y 2 = b and y S exp(f)). 

The y can be computed by replacing the diagonal blocks of b, when we choose a suitable 
f) . We show examples of selections of f) in Section 13.21 

4. Compute to = pyp^ . Here, to € exp(m) because 0(to) = 0(p)0(y)0(p^) = py^p^ = rw . 

5. Compute k — grrv . Then, k always satisfies k £ exp(t) because (to 2 )t = g^0(g) and 

Q(gm if ) — 0(g)Q(m^) = g(m 2 )^m = gm' . 

Steps 2-4 provide a method for computing the square root of a matrix to find to from to 2 . 
After these procedures, we obtain 

g = kpyp^ = kypK (7) 

Here, 0(fc) = k, Q(y) = y\ and Q(p^) = p', so that the decomposition follows the G — KAK 
theorem. All matrices that appear through the algorithm can be described using g and 0. 
This enables us to show that algorithm can automatically reproduce the well-known QFT 
circuit (See, Section |4~2^) . 

3.2 Methods for performing steps 2 and 3 based on KGD 

Above, we provided a method for computing G = KAK decomposition. We did not provide a 
concrete method for computing p, b, and y in steps 2 and 3 in our algorithm. To show examples 
of such concrete methods, we fixed and f). Here, we treat examples that compute KGD. We 
choose as in (0 and show two methods that compute p, 6, y for particular selections of f); 
i.e., f)i = span {\j){j\®ia x \j = 0, • • • ,2™~ 1 -1} and?) 2 = spa,n{\j)(j\®i(<j x ®cr x +o-y(E)o- y )\j = 
0, • • • , 2™~ 2 — 1}. We should note that the second selection is not always possible because it 
demands that all eigenvalues appearing in to 2 in step 1 should be duplicated twice. 

3.2.1 Example 1 

Let be as in |J?jJ|, and f)i as span <8> ia x \ j = 0,- • • , 2"~ 1 — 1}. Then, we compute p, 

b, and y in steps 2 and 3 as follows: 
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(i) Compute eigenvalue decomposition of m 2 . 

Let m 2 = pdjfl be eigenvalue decomposition and \x\ , ■ ■ ■ , /ijv be the columns of p, where 
N = 2™. Then, d is a diagonal matrix in which diagonals have eigenvalues of m 2 and 
all /ij's are eigenvectors of m 2 and mutually orthogonal. 

(ii) Normalize all fij's by 

+ °"nz)M llMi - ^nzMj II 

(iii) For all za,'s that are associated with imaginary eigenvalues, 

(a) let W\ , W 2 , and W3 be sets of vectors such that 

Wi = {uj I a nz Vj = Vj}, W 2 = {fj I <J nz Vj = -Uj}, W 3 = {vj I o nz Uj ^ ±Uj}. 

(b) For each w € W3, compute 

+ W + VnzW _ W (7„ z ffi 
^ = 71 17/ j ^ = 77 m • ( y ) 

\\w + a nz w\\ \\w-a nz w\\ 

Here, || • || denotes the length of a vector. Then, 

• if all elements in Wi and u + are linearly independent, then W\ = W\ U {^ + }; 

• if all elements in W 2 and u~ are linearly independent, then W 2 = W 2 U 

(iv) Repeat steps (a) and (b) for all jUj's that are associated with positive real eigenvalues. 

(v) Repeat steps (a) and (b) for all fXj's that are associated with negative real eigenvalues. 

(vi) Let p = (vi,V2, ■ ■■ , Vff), where v 2 j-i € W\ and v 2 j € W 2 , for j = 1, • ■ • , AT/2. 

The computation procedure follows from Appendix A. Since <r n z v 2j-i — v 2j~i and 
CnzV 2 j = —V2j, we can easily check that Q(p) = p. Then, b = p^m 2 p satisfy 0(6) = 
and b e exp®, where b = Y*l~o <8> i^PO) (0 < Cj < *"). 

(vii) Compute y by replacing all R X {2Q) appears in b with R X {Q). 

Since Q(R X (2Q)) = i4(Xi), #*(Ci) = #x(2(j), then y satisfies y 2 = 6 and y e exp(f)i). 

In step (vi), one may notice that, when R x (ir) appears in b, then we can use Ry(ir/2) instead 
of R x (ir/2) as a replacement rule. 

Since k and p' E are elements of exp(t) that has the specific pattern as shown in Fig. ^ 
they have the following decomposition: 

*=5i 0) ®|0)(0|+^ 1) ®|l)<l| > P f =5f ®|0)(0|+ 5 «®|1)(1|, (10) 

where g[ j \g ( 2 3) £ SU(2 n ~ 1 ) for k,p^ G SU(2 n ) (j = or 1). Here, and gf ] are composed 
of nonzero elements (black squares in Fig. ^) of odd rows, and g[^ and g^ are composed of 
nonzero elements of even rows. 



Yumi Nakajima, Yasuhito Kawano, and Hiroshi Sekigawa 7 



x 



3 



y 



Fig. 2. Image of a decomposition when we choose as in ||§J and a Cartan subalgebra hi as 
span ® i&x | J = 0> ' " > a 71 " 1 - 1}: 

The matrices uses the same notation as in Fig.^to represent the properties. In the quantum circuit, 



the symbol of the control qubit represents the uniformly controlled rotations [5,16]. 



,0-1) 



G SU(4) 



(I, j = 1 or 2) are applied selectively; that is, and g^ are applied when the third qubit is 
whereas g^ and g^ are applied when it is |1). 



Fig.0shows the image of a decomposition in J5J for a three-qubit system; that is, we choose 
as in © and a Cartan subalgebra as f)i. Applying the decomposition in Q recursively 
to elements ffj ,33 e S'J7(2 n_1 ) (j=0 or 1), we obtain a sequence of uniformly controlled 
rotations like in Fig. 13 in Ref. [5], except that R x is used instead of R y in our case. The 
full decomposition of these uniformly controlled rotations into elementary gates has been 
provided by Mottonen et al. [5, 16]. Also, if we change the order of qubits and apply the 
quantum Multiplexor decomposition to k and jr in (J2J), the produced circuit is the same as 
that in Fig. 2 in Ref. [6], except that rotation R y is used instead of R x . Therefore, the number 
of elementary gates needed to compose g £ SU(2 n ) in our method is 0(4"), which is the same 
as in Refs. [4, 6]. 

3.2.2 Example 2 

Here, we show another example of methods for computing p, b, and y in Section 13. II Here, 
we choose f) 2 = span <8 i(cr x <£> a x + <r y ® a y ) \ j — 0, • • • , 2™ _1 — 1}. The Cartan 

involution is the same as in Q . The decomposition of this type is chosen as an example of 
a decomposition of the QFT. In this case, we only change steps (vi) and (vii) in Section 13.2.11 
as follows: 

(vi) Let p = (vi,v 2 , • • • , v N ), for j = 1,- • • , A/4, 

• ^4j-3 G Wi and it is associated with positive eigenvalues, 

• l> ij~2 £ W2 and it is associated with negative eigenvalues, 

• Vij-i € W\ and it is associated with negative eigenvalues, 

• Vij € W2 and it is associated with positive eigenvalues. 

Then, p also satisfies Q(p) — p, and b — p^m 2 p is a block-diagonal matrix, in which the 
diagonals are constructed from the 4x4 matrix 



(1 0^ 

cos(2£j) isin(2(j) 

isin(2Cj) cos(2£,-) 

\0 1/ 



eXp^Cj^CTx ® + a y <g> CT y )). 

Here, the middle part of the above matrix is R X (2Q). 
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(vii) Compute y by replacing R X (2Q) in b with R X (Q). 

4 Decompositions of the QFT using G = KAK decomposition 

In this section, we show that we can automatically reproduce the well known QFT circuit 
using our method. All the matrices that appear through our algorithm can be described 
using the input matrix g and O. In contrast, it is difficult to describe all the matrices that 
appear through the CSD-based algorithm because, as shown in JJJ, the input matrix g has to 
be divided into four square matrices and SVD has to be applied to each partitioned matrix. 
Furthermore, we have to choose a suitable decomposition for each partitioned matrix <?n, g\2, 
521, and 322 such that Eq. JJJ holds, which makes it difficult to formulate U , E, and V in 0J. 

Fortunately, the n-qubit QFT is a very special matrix that has the following property: If we 
permute the order of qubits, then each partitioned matrix can be described using (n— l)-qubit 
QFT. Using the feature, we provide a decomposition of the QFT by CSD-based algorithm. 
This is shown in Section [Ql 

4-1 Notation 

The QFT on n qubits, F n , is a 2" x 2™ matrix such that 



F„ = [J^^-m-^y^ where Wn = «p(^V (11 l 



We define Q n as a 2™ x 2™ permutation matrix: Q n = Xn—i ' ' ' X2X1 > where Xj is the SWAP 
gate applied to the j-th and the k-th qubits. Let H\ = H ® then il 111 is written as 

F n = ^fp"" 1 "o" 1 ^ -1 W = H!D n (I ® Fn-JQn, (12) 

V2 \fn-l — "n-l-Tn-l/ 

where 

\ _ ,. H 2 n-l_ u 

D n =\ q Q J' Mn-i = diag{l,u n ,--- ,u n ). 

This notation follows from Section 4.6.4 in Ref. [17]. 

4-2 Decomposition of the QFT by our method 

Following Section ITT! we compute a decomposition as follows: 

1. Compute m 2 = 9(F%)F n . 

Let S be (I ® F n _i)Q n . Since 6(i?i) = Hi, Q(D n ) = D n , and ai z S = Sa nz , we have 
m = v\z&nzS. All column vectors of S' are then eigenvectors of m 2 because a\ z a nz 
is a diagonal matrix in which diagonal elements are eigenvalues of m 2 . 

2. Decompose m 2 = pbp 1 such that O(p) = p and 0(6) = b< . 

We define p = S^Q n (This selection was done so that p satisfies Q(p) = p and follows 
Section 1X7^ . then 6 = p^m 2 p = Qta lz a nz Q n = /®("" 2 ) ® diag\l, -1, -1, 1). 
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3. Choose y such that 0(y) = y'. Following the step (vi) in Section 13.2.21 we have 

= /®(™- 2 ) (gi exp(n(a x <g> a x + a y ® cr„)/4) 



y = J®(ti-2) 



/l 0\ 
t 
i 
\0 1/ 

This is obtained by replacing each 4x4 diagonal block of v, b = diag(l, —1, — 1, 1), 
with y = exp(Tr(a x ® a x + <r y ® a y )/i). We can easily check that 8(y) = y^ and y 2 = 6, 
because 9(y) = and y 2 = &. 

4. Compute m = pyp^ . 

5. Compute fc = gm* = F n rrv . 

Then, k = kp = (F n py'p')p = F n py', so we have the following decomposition: 

F n = kyp^ = (HiD n Q n y ] )y(Q\ l S) = H x D n S = H t D n {I ® F n ^)Q n . (13) 

We apply a similar decomposition to Fj, for j = n — l,n — 2, •••,2. Next, we show a 
decomposition of D n . D n is controlled-f2j (where j = n — 1, n — 2, • • • , 2), so it suffices 
to consider the decomposition of Qj. Since Qj 6 cxp(t) (it follows from G(fij-) = Qj-), we 
apply the decomposition in l|l(J|) to fl,-. Consider £1$ as an example, then we have O3 = 
gf ] ® |0> <0| +5^ <g> |1)(1|, where ff J 0) = diag (1, w 2 , w 4 , w 6 ) and 5 ^ = rfwg (w, oj 3 , uj 5 , uj 7 ) = 
u> ■ diag (1,uj 2 ,uj 4 ,uj 6 ). Then, we have = diag (1, uj 2 , uj 4 , uj 6 ) (8 e?iag (1,cj). Similarly, since 
diag(l, uj 2 , uj 4 , uj 6 ) is also an element of exp(t), it is decomposed into diag(l, w )®dioflf(l, w 2 ). 
Therefore, fl n _i is composed of n — 1 single-qubit rotations as follows: 

0„_r = (J ^-x) w fO '" (o -J' (U) 

The circuit obtained from the above decomposition is shown in Fig. |3| We apply a similar 
decomposition to Fj (for j = n — 1, n — 2, • • • ,2). Finally, we have a full QFT decomposition 
composed of n Hadamard gates, \n{n — 1) controlled-rotations, and |_§J SWAP gates. Here, 
SWAP gates that appeared in a sequence of permutations Q n Qn-i ■ • • Q2 were optimized. It 
is known that a controlled-rotations can be implemented by three single-qubit rotation and 
two CNOTs, so that the number of elementary gates that appear in Fig- El is 0(n 2 ). 

4.3 Decomposition of the QFT by CSD 

If we apply to the input matrix F n beforehand, then we have the following decomposition: 
F n Ql = UEV 

Sj®n-1 o \ l /J®(n-1) _J®(n-l)\ Q 

r*»-l) ■ ^J®(n-1) j®(n-l) ) ■ y -Q n _iF n _ X/ 

(i?^®7®("- 1 ))(^®/^"- 1 ))£> n (7®F n _ 1 ) = HiS. (16) 



,(15) 



We can easily check that (|15H satisfies the definition of CSD, where all 0/s appear in c and 
s in (IJ are 7r/4. Although we have to apply Qj beforehand when we apply CSD recursively 
to Fj (j = 1, • • • , n — 1), we obtain the well-known QFT circuit. Using the feature, Tucci 
reproduced the well-known QFT circuit using a CSD-based algorithm [2,14]. 
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— R,— 



R" 



H 



■-X-X x-x- 



■ R"~ 



Qn 



D„ 



Fig. 3. Decomposition of the QFT by our algorithm. Here, R J n = diag (l,u^), i.e., a single-qubit 
rotation, and Q n is composed of n — 1 SWAP gates. 



4-4 Decomposition of the QFT by the quantum Shannon decomposition (QSD) 

The QSD [6] is a method that combines CSD and the quantum multiplexor decomposition. 
First, we compute CSD for the input matrix and then apply the quantum multiplexor decom- 
position (cf. [6], Theorem 12), 



ui 0\ _ fv 0\ fd 0\ fw 



(17) 



to {/ and V in |T]l. Here, ui?4 = vd 2 v^ and u> = dv^U2- In the QFT, we apply the decompo- 
sition in I jl 7(1 to y because t7 is an identity matrix in 1(15(1 . Then, we have w = J®^" 1 ) and 

d = \ f— O^Zi) and to = £l„_iF„_i in (TP?! . Here, ^/ — fi n _i is a 2"~ 1 x 2"~ 1 diagonal 

matrix whose (j,j)-th component is ioj^Tj, (j = 1,- • • , 2™~ 1 ). Therefore, the decomposition 
of QFT by the QSD is as follows: 

F n Qn = exp(«tr y <g> S 2 ) exp(-icr 2 <g> 5 3 )(7 ® w 4 ), (18) 



where U4 = f^n-i-^n— 1 an( i ^2 and S3 are 2™~ 1 x 2 n_1 diagonal matrices. Each element of 
62 is 7r/2™ _1 and each (j, j)-th element of $3 is — j7r/2 n+1 . Furthermore, y/Q n -i is composed 
of n — 1 single-qubit rotations as follows: 

Here, ((18|) is also equal to the well-known QFT decomposition l(12|) after optimization as shown 
in Fig. To obtain the decomposition in ((18(1 . note that we have to apply Qj (j = 0, • • • , n) 
beforehand for Fj. 

5 Conclusion 

We introduced a new algorithm for computing any type-AIII KAK decomposition according 
to the given global Cartan involution 0. Recursively performing the decomposition leads us to 
a quantum circuit composed of uniformly controlled rotations. The algorithm can derive any 
matrix decomposition corresponding to the type-AIII KAK decomposition, and it contains 
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D 



D®D 



exp(7 ® v 4 ) ' exp(i c ® S 2 ) 



it* 



Fig. 4. Decomposition of the QFT by the QSD. Here, we move Q n to the right-hand-side by 
inverting it. The dark-gray block in the middle circuit can be simplified. We show the simplified 
block by the same color in the rightmost circuit. Here, D = ■J — fi n _i. 



CSD and KGD as its special cases. This is because our algorithm contains arbitrariness in 
selecting the Cartan subalgebra t) and a square root matrix m for the given Cartan involution 
0, where m is a matrix derived from a global Cartan decomposition g — km. We also showed 
two methods for computing a square root matrix. 

Although the correctness of our algorithm depends on Lie group theory, the main methods 
involved are eigenvalue decomposition and a simple replacement rule. Thus, we can compute 
a decomposition without knowledge of Lie group theory. 

As an example, we showed that our method automatically reproduces the well-known 
QFT circuit for arbitrary n-qubits. When using CSD-based algorithms, we have to apply 
some permutations beforehand in order to reproduce the circuit. The same technique can 
not always be used to describe the canonical form of the decomposition for a given matrix 
because matrices do not always have a convenient form like QFT. Our algorithm might be 
useful in showing the effectiveness of G = KAK matrix decompositions for other particular 
input matrices, because all matrices appearing through our algorithm can be described using 
a given input matrix g and a given Cartan involution 0. 
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Appendix A Properties of eigenvalues of m 2 = Q(G')G 

In Section IS. 2. II we show that 6 is a block-diagonal matrix, each block of which is R x (2(), 
where f) = span{|j)(j| (gi ia x \ j = 0, • • • , 2 n ~ 1 — 1} C m. This appendix provide two lemmas 
and proofs for the eigenvalues of m 2 — 9(G^)G, where G € SU(2 n ) to prove that. Let 
X = m 2 = 9(G t )G. 

Lemma A.l (Properties of complex eigenvalues of X) 

1. The number of the complex eigenvalues of X , with multiplicity counted, is even. They 
are of the form at\, ~ct\, ct2, 02, ■ ■ ■ , on, en, repeated with multiplicity. 



2. There exists unit length vectors u\, . . . , U2t that are mutually orthogonal and satisfy 
o~nzU2j-i = u 2j-i arid <r nz U2j = —U2j- Furthermore, there exists £i, Ct £ M 

satisfying Xu2j-i = cos(Q)u2j-i — ism.(Q)u2j and Xu 2 j = — isin(£j)u2j_i+cos(£j) u 2j- 



Proof. Let ai, . . . , «2» be the eigenvalues of X, repeated with multiplicity. Then, except 
for the order, the elements of the sequence en, . . . , tt^ are equal to ai, . . . , a2>» because 
Q(X) = X\ Therefore, the number of the complex eigenvalues, counted with multiplicity, is 
even, and we denote it as It. We can write the complex eigenvalues as a.\, a±, . . . ,a.t, at- 
That is, the first statement holds. 

Let u be an eigenvector of X corresponding to the eigenvalue a, then a nz u is an eigenvector 
of X corresponding to the eigenvalue a because X(a nz u) = a nz X'u = ctj(a nz u). Now, 
let W be the eigenspace corresponding to a complex eigenvalue a and /3i, . . . , f3 r be an 
orthonormal basis of W. Then, let a nz Pi , . . . , <J nz (3 r be an orthonormal basis of the eigenspace 
corresponding to the eigenvalue a. Therefore, U2j-i = Pj + &nzf3j and u%j — fij — o nz (3j for 
j = !,■■■ , r, span the eigenspace of X corresponding to the eigenvalues a and a. Furthermore, 
they are eigenvectors of o~ nz because o~ nz U2j-i — u-2j-i and o~ nz U2j = —U2j- On the other 
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hand, the following relations hold: 

Xu 2 j-i = af3j +aa nz (3j 
Xu 2 j = af3j - aa nz f3j 
Therefore, put ( = — arg(a). Then cos(C) = Re(a) and sin(£) = — lm(a). That is, 

I«2j-i = cos(C)u 2 j-i - isin(C)u 2j , Xu 2j = -i sm(C)u 2 j-i + cos(()u 2j . 
Similar arguments hold for the other complex eigenvalues. □ 

Lemma A. 2 (Properties of real eigenvalues of X) Real eigenvalues of X are ±1. 

1. The multiplicity of the eigenvalue 1 is even. There exists an orthonormal basis u 2 t+j 
(j = 1, ••• ,2/i) of the eigenspace corresponding to the eigenvalue 1 that satisfies the 
o~nzU 2 t+ 2 j-i — u 2t + 2 j-i and a nz u 2 t+ 2 j = —u 2t + 2 j- 

2. The multiplicity of the eigenvalue —1 is even. There exists an orthonormal basis u 2t + 2fl +j 
(j = 1, •• ■ , 1v) of the eigenspace corresponding to the eigenvalue —1 that satisfies 
o~nzU 2 t+ 2fl + 2 j-i = u 2t + 2 i_ l + 2 j-i and a nz u 2 t+ 2yi + 2 j = — 1i2t+2ju+2i • 

Proof. Lemma [A. II implies that the product of all the complex eigenvalues is 1, and thus 
the product of all the real eigenvalues is 1. Therefore, the real eigenvalues of X are 1 or 

— 1, thus both of the multiplicities of the eigenvalues 1 and —1 are even. Let Wi and W 2 
be the eigenspaces of X corresponding to the eigenvalues of 1 and —1, respectively. Then 
(T nz Wi C Wi and a nz W 2 C W 2 hold because Q{X) = X*. Put W = W x © W 2 . Then WMs 
the direct sum of the eigenspaces for the complex eigenvalues of X. The trace of o~ nz \w 
is because Lemma IA.1I implies that the trace of o~ nz \ w ± is 0, and this implies that the 
multiplicities of the eigenvalues of 1 and —1 of cr nz \w are equal. We write W\ = W\\ © W\ 2 
and W 2 = W 2 i © W 22 , where W\\ and W21 are the eigenspaces of a nz corresponding to the 
eigenvalue 1, and Wi 2 and W 22 are the eigenspaces of a nz corresponding to the eigenvalue 

— 1. Let dim be dij. Then we have dn + d 2 \ — d\ 2 + d 22 . 

We can make similar arguments forX' = 6(G)G t . We write W[ = W[ x © W{ 2 and 
W 2 = W 21 © W 22 , where W xl and W 21 are the eigenspaces of a nz corresponding to the 
eigenvalue 1, and W[ 2 and W 7 ^ are the eigenspaces of o~ nz corresponding to the eigenvalue 
-1. Let dim W-j be d' i} . Then we have d' n + d' 21 = d' 12 + d 22 . 

From Lemma 1a. 31 below, we have dn — d'n, d\ 2 = rf' 12 , d 2 \ = d' 22 , and c?22 = d' 21 . Then, 
we have dn — d\ 2 and d 2 \ — d 22 . Therefore, the statements of the lemma hold. □ 

Lemma A.3 Let X = e(G t )G and X' = 6(G)G t . 

1. If Xu = u and o~ nz u = u, then X'(Gu) — Gu and o~ nz (Gu) = Gu. 

2. If Xu = u and a nz u — —u, then X'{Gu) — Gu and a nz (Gu) = —Gu. 

3. If Xu = —u and o~ nz u = u, then X' (Gu) = —Gu and o~ nz (Gu) = —Gu. 



Ot © OL OL — Ct 

— 2 — U2j-i-\ 2 — U2 i = Re(a)ti2j_i +ilm(a)ti2j, 

OL — ~C£ Ct © ~OL 

— 7. — "2j-i H 7y — u 2j = ilm.(a)u 2 j_i + Re(a)u 2j -. 
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4- If Xu = —u and a nz u = —u, then X'(Gu) = —Gu and <r nz (Gu) = Gu. 

Proof. First, we prove the statements for the eigenvalues and eigenvectors of X' . Since 
G^X'G = @{X) holds, we have (X'(Gu),Gu) = {G^X'Gu,u) = {Q{X)u,u) = {X^u,u) = 
(u,Xu). Thus, the equation Xu — eu, where e = ±1, implies (X'(Gu),Gu) = e(u,u). On 
the other hand, the Cauchy-Schwarz inequality implies \(X'(Gu),Gu)\ < \\X'(Gu)\\ ■ ||Gw||. 
Since the right-hand side is equal to ||Gw|| 2 = ||m|| 2 = (it, u) = \(X'(Gu),Gu)\, the equality 
\{X'{Gu),Gu)\ = \\X'{Gu)\\ ■ \\Gu\\ holds. Therefore, we have X'(Gu) = a{Gu) for some 
a e C; that is, Gu is an eigenvector of X' corresponding to the eigenvalue a. The equality 
a = e follows e(u, u) = (X'(Gu),Gu) — (aGu, Gu) = a(Gu, Gu) — a(u, u). 

To prove the statements for the eigenvalues and eigenvectors of a nz , we use the following 
equations: 

(cr nz (Gu),Gu) = (G^a nz Gu,u) = (<j nz Xu,u) = (Xu,<r nz u). 

The equations Xu = eu and a nz u = e'u, where e, e' = ±1, imply (a nz (Gu),Gu) = ee'(u,u). 
From similar arguments for the eigenvalues and eigenvectors of X', we have <r nz (Gu) = 
ee'(Gu). □ 



